Extreme Outages due to Polarization Mode Dispersion 
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We investigate the dependence of the bit-error-rate (BER) caused by amplifier noise in a linear 
optical fiber telecommunication system on the fiber birefringence. We show that the probability 
distribution function (PDF) of BER obtained by averaging over many realizations of birefringent 
disorder has an extended tail corresponding to anomalously large values of BER. We specifically 
discuss the dependence of the tail on such details of the pulse detection at the fiber output as 
"setting the clock" and filtering procedures. 
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Transmission errors in modern optical telecommunica- 
tion systems are caused by various impairments (limiting 
factors). In systems with the transmisson rate 40 Gb/s 
or higher, polarization mode dispersion (PMD) is one 
of the major impairments. PMD leads to splitting and 
broadening an initially compact pulse [J, 0, La, lil The ef- 
fect is usually characterized by the so-called PMD vector 
that determines the leading PMD-related pulse distortion 
[M El • It is also recognized that the polarization vec- 
tor does not provide a complete description of the PMD 
phenomenon and some proposals aiming to account for 
"higher-order" PMD effects have been recently discussed 
H, El El El- Birefringent disorder is frozen (i.e. it does 
not vary at least on the time scales corresponding to 
the optical signal propagation). Optical noise originat- 
ing from amplified spontaneous emission constitutes an 
impairment of a different nature: The amplifier noise is 
short correlated on the time scale of the signal width. In 
this letter we discuss the joint effect of the amplifier noise 
and birefringent disorder on the BER. Our main goal is 
estimating the probability of special rare configurations 
of the fiber birefringence that produce an anomalously 
large values of BER, and thus determine the information 
transmission reliability. Evaluation of the signal BER 
due to the amplifier noise for a given realization of bire- 
fringent disorder is the first step of our theoretical anal- 
ysis. Second, we intend to study the PDF (normalized 
histogram) of BER, where the statistics is collected over 
different fibers or over the states of a given fiber at dif- 
ferent times, and focus on the probability of anomalously 
large BER. We analyze the basic (no compensation) sit- 
uation and compare it with the case of the simplest com- 
pensation scheme known as "setting the clock". More 
sophisticated compensation strategies will be discussed 
elsewhere. 



The envelope of the optical field propagating in a given 
channel in the linear regime (i.e. at relatively low optical 
power), which is subject to PMD distortion and amplifier 



noise, satisfies the following equation 0, El 113 

d z V - iA{z)V - m{z)d t & - id(z)d%9 = £(z, t). (1) 

Here z, t, £, and d are the position along the fiber, 
retarded time, the amplifier noise, and the chromatic 
dispersion, respectively. The envelope \P is a two- 
component complex field, the two components represent- 
ing two states of the optical signal polarization. The 
birefringent disorder is characterized by two random 2x2 
traceless matrix fields related to the zero-, A, and first-, 
to, orders in the frequency expansion with respect to the 
deviation from the carrier frequency lu . Birefringence 
that affects the light polarization is practically frozen (t- 
independent) on all propagation-related time scales. The 
matrix A can be excluded from the consideration by the 
transformation \P — > V\P, £ — > V£ and rh — > VmV^ 1 . 
Here, the unitary matrix V(z) = Tcxp[i j Q dz'A(z')] is 
the ordered exponential defined as a formal solution of 
the equation d z V = iAV with V(0) = 1. Hereafter, we 
will always use the renormalized quantities. We further 
represent the solution of Eq. (JIJ as \P — tp + <fi where, 

<p = W(z)¥ (t), 0= [dz'W{z)W- 1 (z , )Z(z',t), (2) 



W(z) = cxp 



i / dz'd(z')d; 
I Jo 



Texp 



dz'm(z')dt 



, (3) 



and \Po(t) stands for the initial pulse shape. 

We consider a situation when the pulse propagation 
distance substantially exceeds the inter-amplifier separa- 
tion (the system consists of a large number of spans). 
Our approach allows to treat discrete (erbium) and dis- 
tributed (Raman) amplification schemes within the same 
framework. The additive noise, £ generated by optical 
amplifiers is zero in average. The statistics of £ is Gaus- 
sian with spectral properties determined solely by the 
steady state features of amplifiers (gain and noise figure) 
The noise correlation time is much shorter than 
the pulse temporal width, and therefore £ can be treated 



2 



as ^-correlated in time. Eqs. (121.' il) imply that the noise 
contribution to the output signal <p is a zero mean Gaus- 
sian field characterized by the following pair correlation 
function 

(Mz,h)tf 3 (z,t 2 )) = D^ZS^Sih -t 2 ). (4) 

and, therefore, is statistically independent of both d{z) 
and m(z). Here, Z is the total system length, the prod- 
uct D^Z being the amplified spontaneous emission (ASE) 
spectral density of the line. The coefficient is intro- 
duced into Eq. EJ to reveal the linear growth of the ASE 
factor with Z 15]. 

The matrix of birefringence m can be parameterized 
by a three-component real field hj , rh = "^2 hj&j, with Uj 
being a set of three Pauli matrices. The field h is zero 
in average and short-correlated in z. The above transfor- 
mation rh — » VrhV' 1 guarantees the statistics of hj to 
be isotropic. Since h enters the observables described by 
Eqs. l]2l.'ill in an integral form the central limit theorem 
(see, e.g., implies that the field hj can be treated as 
a Gaussian field with 



(hi(zi)hj(z 2 )) 



D m SijS(zi 



*2 



(5) 



and the average in Eq. 10) is taken over the birefringent 
disorder realizations (corresponding to different fibers or 
states of a single fiber at different times). In case of 
weak birefringent disorder the integral H = J Q Z dz h(z) 
represents the PMD vector. Thus, D m = k 2 /12, where k 
is the so-called PMD coefficient. 

We consider the return-to-zero (RZ) modulation for- 
mat when the pulses are well separated in t. The signal 
detection at the line output, z = Z, corresponds to mea- 
suring the output pulse intensity, /, 



1=1 dtG{t)\K<p(Z,t) + K<t>(Z,t)\ 2 , 



(6) 



where G(t) is a convolution of the electrical (current) fil- 
ter function with the sampling window function. The 
linear operator JC in Eq. © stands for an optical filter 
and a variety of engineering "tricks" applied to the out- 
put signal, \P(Z,t). Out of the variety of the "tricks" 
we will discuss here only those correspondent to optical 
filtering and "setting the clock" compensation. The lat- 
ter can be formalized as, K, c i^ — \P(t — t c i), where t c i is 
an optimal time delay. Ideally, I takes two distinct val- 
ues corresponding to the bits "0" and "1", respectively. 
However, the impairments enforce deviations of I from 
the ideal values. The output signal is detected by in- 
troducing a threshold (decision level), Id, and declaring 
that the signal codes "1" if I > Id and "0" otherwise. 
Sometimes the information is lost, i.e. an initial "1" is 
detected as "0" at the output or vise versa. The BER is 
the rate of such events which is extracted from measure- 
ment of many pulses coming through a fiber with a given 



realization of the PMD disorder, hj(z). For successful 
system performance the BER should be extremely small, 
i.e. typically both impairments can cause only a small 
distortion of a pulse or, stated differently, the optical 
signal-to-noise ratio (OSNR) and the ratio of the squared 
pulse width to the mean squared value of the PMD vec- 
tor are both large. OSRN can be estimated as Ia/(D^Z) 
where 7 = J dt |<?o(t)| 2 is the initial pulse intensity and 
the integration goes over a single slot populated by an 
ideal (initial) pulse, coding "I". 

Based on Eq. © one concludes that the input "0" is 
converted into the output "1" primarily due to the noise- 
induced contribution cf> and, therefore, the probability of 
such event is insensitive to the PMD disorder in accor- 
dance with Eq. (@J. Therefore, anomalously large values 
of BER originate solely from the "I — * 0" transitions. 
Let B be the probability of such an event. Since OSNR 
is large, B can be estimated using Eqs. 1 IKifl as the prob- 
ability of an optimal fluctuation of <fi leading to / = Id- 
Then one concludes that the product D^ZlnB depends 
on the disorder, the chromatic dispersion coefficient and 
the measurement procedure, while being insensitive to 
the noise characteristics. Since OSNR is large, even a 
weak PMD disorder could produce a large increase in 
the value of B. This fact allows a perturbative evalua- 
tion of the independence on the disorder, starting with 
an expansion of the ordered exponential © in powers 
of h. If no compensation is applied the linear term is 
prevailing. It is convenient to introduce a dimensionless 
coefficient /ii in accordance with (D^Z/Io)\n(B/Bo) = 
HiHs/b + 0(H 2 ), where the initial pulse \Po is assumed 
to be linearly polarized, \Pq (b is the signal width), Bq 



being a typical value of B that corresponds to hj 



0, 



hi Bo ~ —Iq/(D^Z). "Setting the clock" compensation 
cancels out the linear in H contribution, i.e. /ii — if t c i 
is chosen to be equal to H%. In this case and also when 
the output signal is not chirped (this e.g. corresponds 
to the case when the initial signal is not chirped and 
the integral chromatic dispersion, J Q dz d, is negligible) 
one gets {D t Z/I Q ) In(B/B ) = ^(Hf + H 2 )/b 2 + 0(H 3 ), 
where \x 2 is another dimensionless coefficient. 

Intending to analyze the dependence of the parameters 
Tq = — (D^Z/Iq) In Bq, hi, and [ii on the measurement 
procedure, we present here the results of our calculations 
for a simple model case. We assume the Lorentzian shape 
of the optical filter: K f V = J °° dt' e?q>{-t' / t)& {t-t') / t . 
Then, as it follows from Eq. Q the statistics of the in- 
homogeneous contribution, K<p, is governed by the PDF, 
V: 

1 



lnp(0) 



D^Z 



dt 



(7) 



The large value of OSNR justifies the saddle-point ap- 
proximation for calculating B. The saddle-point equa- 
tion, found by varying Eq. J7J over cj), is 



r 2 8? 



iG(t)] /C</> = uG(t)JC<p, 



(8) 
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FIG. 1: Dependence of To and of /ni ; 2 on T,t (measured in 
the units set by re-scaling the pulse width to unity). 



where it is a parameter to be extracted from the self- 
consistency condition ©. B can be estimated by V(4>o)i 
with 4>q being the solution of Eqs. I|6I8|I for / = Id. 
Next, we assume that G(t) = 1 at |t| < T and it is zero 
otherwise. Then, for a given value of u, the solution 
of Eq. ijHJ) can be found explicitly. The value of the 
parameter u is fixed implicitly by Eq. JSJ). Thus u (and 
then B) can be found perturbatively in hj, i.e. as«~ 
uq + 5u, Su <C uq, where uq is the solution of the system 
(1618(1 at hj — 0. For the Gaussian shape of the initial 
pulse, <?o cx exp[— t 2 /(2b 2 )], and Id being the half of the 
ideal output intensity (corresponding to ^(Z) = 'Z'o), the 
numerically found dependencies of To, (j,i,2 on r and T 
(measured in the units of pulse width b) are shown in 
Fig. [U 

The PDF S(B) of B, (which appears as the result of 
averaging over many birefringent disorder realizations) 
can be found by recalculating the statistics of Hj using 
Eq. JSJ followed by substituting the result into the cor- 
responding expression that relates B to Hj. Our prime 
interest is to describe the PDF tail corresponding to Hj 
substantially exceeding their typical value \/ D m Z , how- 
ever, Hj remaining much smaller than the signal duration 
b. In this range one gets the following estimate for the 
differential probability S(B) dB: 



a) exp 



2D 



W *(* 

Bo 



u 2 r 2 



™ b)^^ (9) 

B ' °) B l+a ' 



where (a) corresponds to the no-compensation case, (b) 
is related to the optimal "setting the clock" case, and 
a = D(b 2 /(2fi2D m Io). Note, that the result correspond- 
ing to case (b) shows a steeper decay compared to case 
(a), which is a natural consequence of the compensation 
procedure applied. 

Summarizing, our major result is the emergence of the 
extremely long tail © in the PDF of BER. Note that 



Eq. © shows a complex "interplay" of noise and disor- 
der. To illustrate this focal point of our analysis consider 
an example of a fiber line with the parameters Tq = 0.06, 
/.ti = 0.06 and /i2 = 0.12, which is also characterized by 
typical bit error probability, Bq = 10~ 12 , corresponding 
to I /[D^Z] ps 460. Let us also assume that the PMD co- 
efficient, k — y/l2D m , is 0.14 ps/^/km, the pulse width, 
b, is 25 ps, and the length of the fiber, Z, is 2, 500 km, 
i.e. D m Z/b 2 « 6 • 10~ 3 . Then we can find a probability 
for B to exceed, say, 10~ 8 , that is to become at least 4 
orders larger than B . Our answers and J^b) give 
for this probability 10~ 4 and 10~ 6 , respectively, which 
essentially exceed any naive Gaussian or exponential es- 
timate of the PDF tail. Also note that even though an 
extensive experimental (laboratory and field trial) of our 
analytical result would be of a great value, some numeri- 
cal results, consistent with Eq. 10 are already available. 
Thus, Fig. 2a of 0] re-plotted in log-log variables shows 
the relation between In S and In B close to the linear one 
given by Eq. 0)). 
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